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We discuss conditions to be put on mean-field-like theories to be able to describe fundamental 
physical phenomena in disordered electron systems. In particular, we investigate options for a consis- 
tent mean-field theory of electron localization and for a reliable description of transport properties. 
We argue that a mean-field theory for the Anderson localization transition must be electron-hole 
symmetric and self-consistent at the two-particle (vertex) level. We show that such a theory with 
local equations can be derived from the asymptotic limit to high spatial dimensions. The weight of 
the diffusion pole, i. e., the number of diffusive states at the Fermi energy, in this mean-field theory 
decreases with the increasing disorder strength and vanishes in the localized phase. Consequences 
of the disclosed behavior for our understanding of vanishing of electron diffusion are discussed. 
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I. INTRODUCTION 

Mean-field theories play an important role in the de- 
scription of thermodynamic systems. They are intended 
and used as a first approximation offering a qualitative 
picture of the physics of the studied phenomena. The 
mean-field concept has developed from its initial intu- 
itive ideas of van der Waals and Weiss through the Lan- 
dau theory of critical phenomena to its present sophisti- 
cation and systematics provided by the limit to infinite- 
dimensional lattice models. At present, a modern mean- 
field theory is no longer a weak-coupling approximate 
treatment neglecting spatial fluctuations. It represents 
a comprehensive theory providing a phase diagram in 
the whole range of the input parameters and simulating 
the exact behavior in specific limiting situations. With- 
out a mean-field theory we are mostly unable to identify 
the relevant fluctuations the mean values of which are 
reflected by thermodynamic (order) parameters. Mean- 
field theory is particularly important for critical phenom- 
ena with divergent correlation functions, where it allows 
us to handle singularities in a consistent and manageable 
way and to select the proper low-temperature phase, at 
least above the lower critical dimension. 

Mean-field theories were primarily developed for col- 
lective phenomena in interacting systems. Nontrivial and 
sometimes not easily understandable effects are, however, 
also induced by randomness. Randomness, in connection 
with interaction or with quantum interference, can cause 
significant and sometimes even unexpected changes in the 
behavior of the system. Since mostly no exact solutions 
are available for disordered systems, a mean-field approx- 
imation has become one of the most powerful tools to 
handle fluctuations in the chemical composition of solids. 

Milestones of a mean-field theory for disordered (non- 
interacting) electron systems were laid at the end of the 
sixties and the beginning of the seventies of the last cen- 
tury. The so-called Coherent Potential Approximation 
(CPA) developed at that time is a self-consistent approx- 
imation describing rather accurately the electronic struc- 



ture and thermodynamic properties of random alloys not 
only at the model level but also in realistic settings. 1 
Later on, the CPA was shown to fit the modern defini- 
tion of the mean-field theory as an exact solution of the 
model system in infinite spatial dimensions'**^ 

The coherent potential approximation is nowadays 
considered as a paragon for mean-field theories of quan- 
tum disordered and interacting systems. Its generalized 
form4 offers one possible interpretation of equations of 
motion in the Dynamical Mean-Field Theory (DMFT)£ 
In spite of the proved reliability of the CPA to produce 
an accurate equilibrium electronic structure of disordered 
systems, it fails in encountering for inter-site quantum 
coherence and backscattering effects. The CPA is es- 
sentially unable to go beyond the semi-classical descrip- 
tion of transport properties qualitatively captured by the 
Boltzmann equation. This inability is due to the fact 
that the CPA does not include vertex corrections to the 
electrical conductivity independently of how strong the 
disorder may be£ It is hence unsuitable for the descrip- 
tion of one of the most prominent features of disordered 
systems: Anderson localization. 

Anderson localization in disordered or amorphous 
solids takes place when there are available electronic 
states at the Fermi surface but no diffusion or charge 
transport at long distances is observed. Possibility of the 
absence of diffusion in impure metals and alloys was pro- 
posed by P. W. Anderson on a simple tight-binding model 
of disordered noninteracting electrons^ Since then, van- 
ishing of diffusion, now called Anderson localization, 
has attracted much attention of both theorists and ex- 
perimentalists In spite of a considerable portion of 
amassed experimental data, disclosed various specific and 
general aspects of the Anderson metal-insulator transi- 
tion, and a number of theoretical and computational ap- 
proaches so far developed we have not yet reached com- 
plete understanding of Anderson localization. Although 
many features of the critical behavior at the Anderson 
localization transition have been disclosed, the position 
of this disorder-driven metal-insulator transition within 
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the standard classification scheme of phase transitions 
with control and order parameters has remained unclear. 
It has been mainly due to the nonexistence of an appro- 
priate mean-field-like theory for this phenomenon. Only 
very recently the present authors demonstrated that a 
mean-field-like solution for the Anderson metal-insulator 
transition can be derived from the asymptotic limit to 
high (but finite) spatial dimensions^ This solution is 
very close to the CPA in the resulting electronic struc- 
ture, it reduces to the CPA in infinite dimensions. It, 
however, differs from the CPA significantly in transport 
properties derived from two-particle functions. In addi- 
tion to the one-electron self-consistence of the CPA, the 
new mean-field theory is endowed with a two-particle self- 
consistence. That is, the two-particle irreducible vertices 
are determined from self-consistent nonlinear equations. 

The mean-field-like theory for the disorder driven van- 
ishing of diffusion of Ref . shows some unexpected fea- 
tures. It contradicts the widely extended dogma that the 
weight of the diffusion pole, i. e., the number of diffus- 
ing particles, does not depend on the disorder strength. 
The weight of the diffusion pole is conserved if all the 
states near the Fermi energy are finite combinations of 
Bloch waves for any configuration of the random poten- 
tial. Or equivalently if a Ward identity between self- 
energy and the irreducible electron-hole vertex holds for 
all transfer energies liLiS We demonstrated in Ref. [l3| 
that once the electron-hole irreducible vertex contains 
the so-called Cooper pole, the number of extended states 
at the Fermi level decreases with increasing the disorder 
strength. Hence, there is no chance for a theory with the 
Cooper pole to fully satisfy the Ward identity between 
the averaged one- and two-electron Green functions and 
to keep the number of diffusive states independent of the 
disorder strength. 

The aim of this paper is to clarify the ambiguities 
connected with the mean-field concept applied to two- 
particle functions in disordered electron systems. We de- 
limit the content and the range of validity of the two 
existing mean-field theories for the Anderson model of 
noninteracting electrons — the CPA and that of Ref. ^3- 
The common ground for both theories is the limit to high 
spatial dimensions. We show that the ambiguity in the 
identification of the mean-field theory for two-particle 
functions results from two different ways how it can be 
derived: either via the generating local functional and 
Ward identities, or via a direct diagrammatic construc- 
tion in high dimensions. The former construction, CPA, 
is suitable only for one-electron spectral and thermody- 
namic properties. Since the CPA lacks the electron- hole 
symmetry and does not correctly reproduce the high- 
dimensional limit of two-particle (vertex) functions, it 
becomes unreliable when applied to the electrical con- 
ductivity and other transport properties. We demon- 
strate that the latter approach, when properly formu- 
lated, can lead to a theory being self-consistent and man- 
ifestly electron-hole symmetric at the one- and as well at 
the two-particle level. The last two conditions are neces- 



sary ingredients for a theory being able to describe the 
Anderson localization transition. 

The layout of the paper is as follows. In Section [H] 
we summarize basic properties of the CPA defined from 
the limit to infinite spatial dimensions. We show how the 
averaged grand potential is derived from the local one- 
particle propagator and the self-energy. The higher-order 
Green functions are then determined via local external 
perturbations. The mean-field theory with a two-particle 
self-consistence is constructed in Section 11111 First, in- 
ability of the CPA to reproduce the proper infinite- 
dimensional limit for two-particle Green functions is 
demonstrated. Then, using the parquet scheme and the 
electron-hole symmetry we derive a self-consistent (non- 
linear) equation for the irreducible two-particle vertex. 
This equation is then solved at the mean- field level, i. e., 
in the leading nontrivial order of the high-dimensional 
limit. The explicit form of the diffusion pole in high spa- 
tial dimensions with its weight dependent on the disorder 
strength is finally obtained. Consequences of our findings 
for understanding of the disorder-driven vanishing of dif- 
fusion and Anderson localization are discussed in the last 
Section HVI 

II. THERMODYNAMIC MEAN-FIELD 
THEORY: ONE-PARTICLE SELF-CONSISTENCE 

A. One-particle functions and generating 
functional 

To construct a comprehensive mean-field theory for 
thermodynamic properties of a random system means to 
find an approximate representation in closed form for the 
grand potential averaged over random configurations 

n(/x) = -i(lnTrexp{-/3^ + /3/xiv})^ (I) 

where /i is the chemical potential and N is the particle 
number operator. We will consider in this paper only a 
noninteracting lattice electron gas scattered on random 
impurities and described by the Anderson tight-binding 
Hamiltonian 

<ij> i 

where Vi is a local, site-independent random potential. 

It has become evident since the introduction of the 
concept of the limit to infinite spatial dimensions in quan- 
tum itinerant systems that a controllable comprehensive 
mean-field theory of itinerant models should be defined 
via this formal limit i 4 In high spatial dimensions the 
diagonal (local) and off-diagonal (nonlocal) elements of 
the one-particle propagator separate from each other. 
The former are of order O(d ), while the latter vanish 
as dr 1 / 2 , where d is the spatial dimension. The full one- 
particle propagator and the self-energy have the following 
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high-dimensional asymptotics 

G = G dia9 [d°} + G off [d- 1/2 } , 
E = E*" 9 ^ ] + E o// [cT 3/2 ] . 



(3a) 
(3b) 



We can classify contributions to the many-body per- 
turbation expansion for the self-energy according to their 
high-dimensional asymptotic contribution and obtain in 
the leading order a local approximation for the irre- 
ducible part of the one-electron propagator. The inter- 
acting part of the thermodynamic potential in infinite 
spatial dimensions is then a functional of only G dla9 and 

In disordered systems the inter-particle interaction is 
replaced by correlations between scatterings on impuri- 
ties. The self-energy is here a coherent potential of an 
effective homogeneous (nonrandom) medium represent- 
ing the effect of impurity scatterings on the motion of 
electrons. Since the scatterings are static, we can find an 
explicit representation of the averaged grand potential in 
infinite spatial dimensions. We can write 



G, E 



F ^Q dia 9 -1 _|_ g<iiag| _ _L ^jy m Qdiag 



Trin 



P 

q(Q)-1 _ j^diag , 



where we denoted 
F 



{x} = --(T,h, 



x-v 



(4a) 



(4b) 



the local " interacting part" of the thermodynamic poten- 
tial, in this case the effect of multiple scatterings^ The 
trace operator Tr extends over the lattice space as well 
as over the Matsubara frequencies. The only nonlocal 
contribution to the generating functional f2 M comes from 

the bare propagator G^ ' . The site-diagonal (local) com- 
plex vectors G dla9 (fi) and E*° 9 (/z) in fermionic Matsub- 
ara frequencies (2n + l)7r//3 are variational parameters, 
the physical values of which are attained at stationarity 
points of the generating functional 0} . 

The defining equation for the local element of the av- 
eraged one-particle propagator is obtained from an equa- 
tion <5£! M /<5E* a5 (/i) = and the self-energy is determined 
from Sfl^/SGf^di) = 0. After straightforward manipu- 
lations the two equations reduce to 



1 + [£„(/*) - Vi\ G„(/x) 



(5a) 



and 



dep{e) 



H - E„(it) - e 

(5b) 

We dropped the superscript diag in the local functions 
and introduced the electronic density of states p(e). Due 
to the static character of the impurity scatterings the 



stationarity equations are diagonal in the Matsubara fre- 
quencies and can be solved for each frequency indepen- 
dently. Inserting the solution for all Matsubara frequen- 
cies to Eqs. (0J we obtain the equilibrium thermody- 
namic potential for noninteracting electrons scattered on 
random impurities. The equilibrium thermodynamics of 
the systems is then determined only by the local irre- 
ducible part of the averaged one-particle resolvent. This 
irreducible part is self-consistently determined from the 
Soven equation J^}. 



B. External sources and two- particle functions 

Thermodynamics of disordered systems is not of much 
interest unless inter-particle interactions are present. But 
even then the averaged thermodynamic potentials de- 
pend on only one-electron functions. One-electron func- 
tions, however, do not contain the complete information 
about the behavior of statistical ensembles, in particu- 
lar of disordered systems. The equilibrium thermody- 
namic potentials do not contain sufficient information 
from which we could derive transport properties of the 
system and its response to weak external electromagnetic 
perturbations. To include the electrical conductivity into 
the mean-field description, the thermodynamic construc- 
tion from the preceding subsection must be extended to 
include averaged two-particle propagators. 

Averaged two-particle propagators in disordered sys- 
tems contain at least two energy arguments (two in non- 
interacting and three in interacting systems). The best 
way to guarantee that one- and two-particle functions 
are approximated consistently within a single approxi- 
mate scheme is to use the Baym-Kadanoff concept of ex- 
ternal sources added to the equilibrium thermodynamic 
potential^ To introduce higher-order Green functions 
with several energies (chemical potentials) into the ther- 
modynamic description we have to replicate the original 
system so as for each energy we have an independent 
replica of the original system, that is, of creation and 
annihilation operators. 

We replicate the creation and annihilation operators 
and introduce external perturbations into the thermody- 
namic description via a generalized grand potential of a 
j/-times replicated system ft v (Ei,E2, . ■ . , E v ; U) with v 
chemical potentials E\, . . . , E v . An external perturba- 
tion U is used to couple different replicas and to break 
the initial replica independence. We then can write 



W{E l ,E 2 ,...,E u] U) 



= -i^lnTrcxp j-0 £ (h%5 13 

+Atf(«))}) (6) 

where we assigned to each replica characterized by energy 
(chemical potential) E{ a separate Hilbert space and de- 
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noted AflW) = J2u U kP^k ] A i] an external perturba- 
tion to be set zero at the end. Thermodynamic potential 
tt u (E 1: E 2 , ...,E V ;U) is a generating functional for aver- 
aged products of Green functions up to the ^th order. In 
practice, we will use linear-response theory with one- and 
two-particle Green functions, i. e., W(E\, E2, ■ ■ ■ , E v \ U) 
is expanded up to U 2 . Therefore it is sufficient to intro- 
duce only two replicas. 

In fact we are interested only in corrections to the 
products of the averaged one-particle propagators ex- 
pressed via vertex functions. The two-particle vertex T 
is defined from the two-particle resolvent G^ in momen- 
tum representation as 

22; q) - G(k, Zl )G(k + q,z 2 ) [<5(k - k') 
+r kk ,(zi,z 2 ;q)G(k',z 1 )G(k' + q,z 2 )] . (7) 

The external disturbance U mixes different replicas 
and propagators in the replicated space are matrices in 
the replica indexes. Since we are interested only in the 
averaged two-particle functions, we can resort to two en- 
ergies and to a two-by-two matrix propagator 



G- 1 (ki, zi, k 2 , z 2 ; U) = G^- 1 + U - S 

'zt-e^-SniU) C/-S 12 (C/) 

U-Z 2 i(U) z 2 -e(k 2 )-£ 22 (£7) 



(8) 



where e(k) is the lattice dispersion relation and the self- 
energy elements Ejj generally depend on both energies 
zi,z 2 . The matrix G represents the averaged resolvent 
that is to be used in the grand potential fl 2 (Ei, E 2 ; U) 
from Eq. 10. It is now a straightforward task to derive a 
matrix Soven equation generalizing Eq. i|5a|) . We obtain 



G(z 1 ,z 2 ;U) 



G- 1 (z 1 ,z 2 ;U) + X(z 1 ,z 2 ;U)-V i 



(9) 



where G(zi,z 2 ;U) = N 2 J2 klk2 G(k 1 ,z 1 ,k 2 ,z 2 ;U) is 
the local element of the matrix one-particle propaga- 
tor. Inversions in Eq. have matrix character in the 
replica space. The diagonal elements of the matrix equa- 
tion © determine the one-particle propagators for ener- 
gies z\ and z 2 . The off-diagonal elements, proportional 
to the perturbation U , determine the local two-particle 
resolvent that is defined as the coefficient at the linear 
term in the expansion of the local matrix propagator 
G(z\,z 2 ]U) in the external perturbation U. The local 
two-particle Green function can be represented with the 
aid of the irreducible vertex (two-particle self-energy) A 
via a Bethe-Salpeter equation. We find from Eq. © that 
the Bcthc-Salpctcr equation in the mean-field approxi- 
mation reduces to an algebraic one 



l{z\,z 2 ) = 



A(zi, z 2 ) 



X( Zl ,z 2 )G(zi)G(z 2 ) 



(10) 



where 7 is the local part of the two-particle vertex T. The 
irreducible vertex A in equilibrium (U = 0) determined 
via Eq. HI U| l obeys an equation 



X(zi,z 2 ) 



5T, u (zi,z 2 ) 



8Gu(zi,z 2 ) 



1 



u=o G(zi)G(z 2 ) 
1 



1 + [£(«!) - Vi] G( Zl ) 1 + [E(z 2 ) - Vi] G(z 2 )/ av 



(11) 



We can easily verify that this equation coincides with the 
CPA solution for the irreducible vertex X(z ll z 2 ). 6 

There is no ambiguity in the mean-field construction 
of local one- and two-particle functions. But a mean- 
field treatment has a physical relevance only if it is able 
to produce nonlocal correlation functions, the long-range 
fluctuations of which may significantly influence the ther- 
modynamic and dynamical behavior. There is not, how- 
ever, a unique way how to generate the two-particle ver- 
tex with non-local contributions within the local mean- 
field approach. The simplest and most straightforward 
way is to use the Bethe-Salpeter equation with the CPA 
irreducible vertex A, Eq. (|10|) . and to replace the local 
propagators with the full nonlocal one-electron propaga- 
tors G(k, z). Such a Bethe-Salpeter equation remains 
algebraic in momentum representation and results in a 
two-particle vertex with only one transfer momentum. 
We obtain 



r ± (z 1 ,z 2 ;q ± ) = 



X(zi,z 2 ) 



1 - \(zi,z 2 )x ± {zi 1 z 2 ;q ± ) 



(12) 



where we denoted the two-particle bubble 

X ± (z 1 ,z 2 ;q) = l^G(k,2 1 )G(q±k,z 2 ) . (13) 



The ambiguity in this definition of the full mean-field 
vertex is in the type of nonlocal multiple scatterings we 
include into the Bethe-Salpeter equation. They are here 
denoted by the superscript ±. The plus sign corresponds 
to multiple scatterings of electron-hole pairs, while the 
minus sign to electron-electron pairs. In case of elastic 
scatterings the electron-hole and electron-electron bub- 
bles produce numerically the same number. However, 
the difference between the two types of pair scatterings 
lies in the respective transfer momentum q*. Using the 
notation for momenta in the two-particle resolvent from 
Eq. J7J we have q + = q and q = q + k + k'. 

This ambiguity in the definition of the mean-field two- 
particle vertex is not usually acknowledged in the lit- 
erature, since the electron-hole scattering channel, rele- 
vant for the electrical conductivity, is preferred and di- 
rectly derived from the Baym-Kadanoff approach. 6 How- 
ever, when the mean-field theory is viewed upon as the 
limit to infinite spatial dimensions, both electron-hole 
and electron-electron multiple scatterings possess the 
same high-dimensional asymptotics^i A priori, neither 
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electron-hole nor electron-electron multiple scatterings 
should be discarded. The appropriate form of the ver- 
tex is then selected by the physical quantities in which it 
appears, such as is the case of the electrical conductivity. 

Incapability of the thermodynamic mean-field the- 
ory to determine uniquely the nonlocal part of the two- 
particle vertex results from the degeneracy of the local 
theory with elastic scatterings only (noninteracting sys- 
tems). Multiple single-site scatterings, the only ones rele- 
vant in the mean-field approach, are unable to distinguish 
between electrons and holes. Only if we include explicitly 
scatterings on distinct lattice sites we are able to distin- 
guish between electrons and holes. Hence, the standard 
thermodynamic mean-field theory of quantum itinerant 
systems uniquely defines only the local two-particle ver- 
tex, while it remains ambiguous in the determination of 
the full nonlocal two-particle vertex. 



III. MEAN-FIELD THEORY FOR VERTEX 
FUNCTIONS: TWO-PARTICLE 
SELF-CONSISTENCE 

A. Nonlocal contributions to the vertex function 

A rather inaccurate way to the momentum-dependent 
two-particle vertex is not the only imperfection of the 
thermodynamic mean-field theory. This theory com- 
pletely fails to account for backscattering effects, ver- 
tex corrections to the electrical conductivity, and the 
Anderson metal-insulator transition. All these effects 
are induced by strong nonlocal quantum coherence and 
spatial correlations reflected in the momentum behav- 
ior of the two-particle (vertex) functions. To capture 
these phenomena we have to go beyond a perturbative 
description and to employ a self-consistent scheme for 
the (irreducible) vertex functions. The best local ap- 
proximation for the irreducible vertex, CPA, is non-self- 
consistent at the two-particle level. We hence have to go 
beyond the CPA and include nonlocal (long-range) con- 
tributions to the vertex function in a non-perturbativc 
manner. Thereby a question arises whether we are able 
to reach a reasonably simple approximation with a two- 
particle self-consistence that could be called a mean-field 
theory. It is clear that such a theory must be momentum 
dependent, but the momentum dependence should be re- 
duced to a necessary minimum. We will demonstrate in 
the next subsections that the desired momentum depen- 
dence can be very effectively reduced by the asymptotic 
limit to high spatial dimensions. 

The best way to construct a mean-ficld-likc approx- 
imation for momentum-dependent functions is to build 
up the theory within a self-consistent expansion around 
the CPA. If we denote the local CPA one-particle propa- 
gator G loc (z) — N^ 1 J2k G loc (k, z), the small parameter 
controlling the expansion around the CPA is a perturbed 
propagator G(k, z) = G(k, z) — G loc (z), where G(k, z) is 
the full one-electron propagator. The CPA propagator 



G loc contains the local self-energy E ioc (z) from Eq. (|5af) . 
while the full one a self-energy S(k, z) that is to be deter- 
mined later from a Dyson equation. We apply the expan- 
sion around the CPA to two-particle functions, where the 
one-particle propagators will be treated as external func- 
tions. It means, that we first disregard the consistence 
between the one- and two-particle functions. This con- 
sistence will be restored later via Ward identities once a 
suitable approximation for the vertex functions has been 
fixed. 

We can classify nonlocal contributions to the two- 
particle vertex by the type of the correlated two-particle 
propagation. We either simultaneously propagate an 
electron and a hole or two electrons (holes). Diagram- 
matically it means that we connect spatially distinct two- 
particle scattering events with antiparallel or parallel 
pairs of one-particle propagators. Multiple scatterings of 
pairs of the same type define a channel of a two-particle 
irreducibility. We call a diagram two-particle irreducible 
if it cannot be split into separate parts by cutting simul- 
taneously either an electron-hole or an electron-electron 
pair of propagators. The two definitions of the two- 
particle irreducibility lead to topologically nonequiva- 
lent irreducible functions and to different Bethe-Salpeter 
equations for the full vertex. In each Bethe-Salpeter 
equation the two-particle functions are interconnected 
via one-particle propagators in a different manner. We 
can generically represent the channel-dependent Bethe- 
Salpeter equations as 

I\k' (z+, z- ; q) = A£ k , (z+ , z_ ; q) 

+ [A"GG0r] tk ,(z + ,z_;q). (14) 

We used the symbol for the channel-dependent mul- 
tiplication of the two-particle functions represented by 
specific momentum convolutions. Here A Q is the irre- 
ducible vertex in the a-channel. 

We will specify the momentum convolutions in the 
generic Bethe-Salpeter equation l|14|l for electron-hole 
and electron-electron multiple scatterings only. There is 
also a third two-particle channel, the so-called vertical 
channel with one-particle self-correlating scatterings^ 
These scatterings are, however, unimportant for the phe- 
nomenon of Anderson localization, since the correspond- 
ing two-particle propagator does not contain the diffusion 
pole. 

Using the notation from Eq. for the momentum 
dependence of the two-particle functions we can write 
explicitly the Bethe-Salpeter equation in the electron- 
hole channel with barred functions as 

I\k' (z+ , z- ; q) = A{&, (z+ , z_ ; q) 

+ ^ E A kk»(^+, z-l q)G + (k")G_(k" + q) 

k" 

x IV k /(z + ,z_;q) . (15a) 
The Bethe-Salpeter equation with the electron-electron 
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multiple scatterings then analogously reads 
IW(z + ,z_;q) =Aw(z+,Z-;q) 

+ ^ E A kk" (*+,*-; q + k' - k")G+ (k") 

k" 

x G_(Q-k")rW(z+,z_;q + k-k") , (15b) 

where wc denoted Q = q+k+k' the transfer momentum 
between the two electrons of the scattered correlated pair. 
In these equations we abbreviated G(k, z±) — * G±(k). 

Equations (|15() constitute the fundamental building 
blocks for the construction of systematic approximations 
for the two-particle vertex. They are analogues of the 
Dyson equation and enable us to determine the full vertex 
from irreducible vertices. We hence can apply the per- 
turbation (diagrammatic) expansion to the two-particle 
irreducible vertices, i. e., two-particle self-energies. As a 
first step toward a mean-field-like theory for these ver- 
tices we have to maximally simplify the momentum de- 
pendence of the vertex functions, but still staying beyond 
the local CPA. This will be achieved by the asymptotic 
limit to high (finite) spatial dimensions. 



The standard nonlocal CPA vertex can be recovered 
from the above expression if we neglect the contribution 
from the electron-electron multiple scatterings, the sec- 
ond term in the brackets on the right-hand side (r.h.s.) 
of Eq. I|18fl . There is, however, no reason for this neglect, 
since both terms within the brackets on the r.h.s. of 
Eq. (|18fl produce the same asymptotic behavior in powers 
of the inverse dimension, 0(d^ 1 ). Their only difference 
is in the momentum dependence. 

The limit to infinite spatial dimensions reduces the 
momentum dependence of two-particle functions but 
does not lead to a mean-field-like approximation for ver- 
tices. There is no self-consistence in the two-particle pa- 
rameters and the local irreducible vertices are determined 
from the CPA. We hence cannot expect that this high- 
dimensional two-particle non-self-consistent asymptotics 
would lead to major deviations from the CPA. The only 
significant change in the vertex function, Eq. (|18f) . with 
respect to the CPA vertex from Eq. (jl2[) are the vertex 
corrections to the electrical conductivity in the form of 
weak localization properly described by multiple electron- 
electron scatterings (maximally crossed diagrams). 



B. Asymptotic limit to high lattice dimensions 

Bcthe-Salpeter equations lj*H)|) use only off-diagonal 
one-particle propagators and hence are suitable for per- 
forming the limit to high spatial dimensions. We use 
the hypercubic lattice that has a straightforward high- 
dimensional limit. The one-electron propagator G has 
the following asymptotics 

S(^)^±t r ^)J {z _^f_ er , (16) 

where E(z) is the CPA (d — oo) self- energy^ 

The irreducible two-particle vertices must collapse 
to local quantities in the limit to infinite dimensions. 
Since the Bcthe-Salpeter equations l|15|) use only the off- 
diagonal propagators vanishing in the limit d = oo, both 
vertices A eh and A ee must coincide with the full local 
two-particle CPA vertex in d = oo. We hence have 

A^(z + ,z-;q) = A^(z + ,Z-;q)=j(z+,Z-) . (17) 

We further denote x(z+,z_;q) = x( z +, z-;q) — G + G_ 
with x(z+,z-;<l) = X + (z+,z-;q), G+ = G loc (z + ), and 
G_ = G (oc (z_). If we take into account only the electron- 
hole and electron-electron multiple scatterings we can 
represent the leading asymptotics of the full two-particle 
vertex in high dimensions as follows 

r kk /(z + ,z_;q) =7(z+,z_) + X(z+,Z-) 
1" l(z+,z-)x(z+,z-;q) 
1 - A(z + ,z_)x(z + ,z_;q) 

l- \(z+,z-)x{z+,z-;Q)\ 



C. Parquet equations and electron- hole symmetry 

To go significantly beyond the CPA predictions for 
transport properties and response functions we have to 
introduce a self-consistence that would extend also to the 
two-particle irreducible vertices. That is, the two-particle 
irreducible vertices A eh and A ee should be determined 
from nonlinear equations. This effect can be achieved by 
introducing the so-called parquet equations. The concept 
of parquet equations is based on the observation that 
two-particle reducible diagrams in one scattering chan- 
nel are irreducible in the other distinguishable scattering 
channels. Parquet equations were introduced in many- 
body theoriesi2i22i2i but recently they were adjusted also 
to disordered systems^ The reducible diagrams from 
one channel can become irreducible in the other channels 
only if different channels are indeed distinguishable or 
nonequivalent. The idea of parquet equations cannot be 
applied to local propagators of noninteracting particles 
with indistinguishable electrons and holes, hence within 
the CPA. It, however, works very efficiently for nonlo- 
cal vertex functions in the Anderson model of disordered 
electrons. 

If we again take into account only the electron-hole 
and the electron-electron scattering channels, wc can 
write the basic parquet equation for the full two-particle 
vertex in high dimensions 

T k k' {z+ , z_ ; q) = A£k' (z+ , z-\ q) 

+ A k e k ,(z + ,z_;q)- 7 (z + ,z_) . (19) 

Equation (|19fl tells us that the full vertex is decomposed 
into irreducible and reducible diagrams in the either scat- 
tering channel and that the reducible diagrams consist 
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of only the irreducible diagrams from the other chan- 
nel from which the completely irreducible vertex, i. e., 
the vertex irreducible in both channels, was subtracted. 
The limit to high lattice dimensions then determines the 
completely irreducible vertex to be the full local CPA 
(d = oo) vertex 7. 

Parquet equation (fH^I can now be used in the Bethe- 
Salpeter equations (|T5)l to exclude the full vertex T from 
them. Thereby we reach a closed set of nonlinear integral 
equations for the irreducible vertices A eh and A ee . This 
set of equations is generally not solvable without further 
approximations. To approximate the resulting parquet 
equations in a systematic and controlled way we again 
utilize the mean- field idea — limit to high dimensions. 
For two-particle functions and parquet equations we have 
to use this limit only in the asymptotic sense so that 
nonlocal fluctuations do not go lost completely. 

One can make an important observation in high spa- 
tial dimensions. The off-diagonal one-particle propa- 
gators G behave in the leading asymptotic order as 
Gaussian random variables with respect to momentum 
summations^ Using representation (|l(j|) we can easily 
prove the following relations 



i^xfa' + q^q' + k^l-G^q-k) , (20a) 
q' 

V + liMq + ^2) = ^ - Q2) , (20b) 



N 



where we used abbreviations Z = t 2 (G 2 + ) (G 2 _) with 
{G±) = N- 1 £ k G±(k) 2 . The equalities in Eq. gD) hold 
only within the leading asymptotic order d — ► 00. The 
functions G± and x form a closed set of Gaussian ran- 
dom variables with respect to momentum convolutions. 
We hence can use Eqs. (l2"Uj) to simplify the parquet equa- 
tions for the irreducible vertices A eh and A ee . 

Before we attempt to resolve the parquet equations 
in high dimensions, we utilize the time- reversal invari- 
ance of the system. It is an important feature of electron 
systems without spin- and orbital-dependent scatterings. 
According to this invariance the physical (measurable) 
results should not depend on the orientation of propa- 
gators. We hence can write for one- and two-particle 
propagators 

G(k,z) = G(-k,z) , (21a) 
r kk /(z + ,z_;q) = r kk /(z + ,z_;-Q) 

= r_ k ,_ k (z+,z_;Q) . (21b) 

In case of the two-particle vertex, the time reversal 
was applied only to one fermion propagator. The 
time reversal leaves the full two-particle vertex invari- 
ant but it transforms the Bethe-Salpeter equation (1 15 all 
to Eq. I|15b|) and vice versa. It means that the irreducible 
vertices transform as follows 



A k e k ,(z + ,z_;q) 



A&(z+,z_;-Q) 



= A' 



l kk' 

eh 



,(z + ,2_;Q) 



(21c) 



The time-reversal (electron-hole) symmetry reduces 
the number of parquet equations to just one nonlinear 
integral equation for a single vertex that we define as 

-Q) 



A kk ,(z+,z 7 ;q) = A-,(z + ,z r ;q) = A&,(z 
The resulting equation for this vertex reads 



+ > • 



Akk'(q) = 7 
+ ^ J2 A kk"(-q k k")G+(k")G_(q + k") 



[A k „ k ,(-q-k'-k") + A k „ k ,(q)- 7 ] . (22) 



Equation (|22|l can now be simplified in high spa- 
tial dimensions by using the Gaussian summation rules, 
Eqs. H20fl . It is clear from these rules that the fermionic 
momenta from different two-particle functions must be 
summed independently in the leading asymptotic order. 
Any correlated momentum summation involving two dif- 
ferent two-particle functions costs a factor l/d. Then 
only the conserved bosonic transfer momenta survive as 
in the case of the high-dimensional vertex from Eq. (|18|l . 
We hence have to sum both sides of Eq. ill'L'l) over in- 
coming and outgoing fermionic momenta k, k' so as to 
extract the high-dimensional limit of the irreducible ver- 
tex A. To reach an equation for the relevant variables we 
introduce 



kk' 



(q) 



Further on we have in the leading order 

jjs £ A kk , (q + k + k' ) = 1 £ A(q) = A . 



(23a) 



(23b) 



kk' 



Since the fermionic momenta from different two-particle 
functions are summed independently in high spatial di- 
mensions, the parquet equation (1221) reduces with the 
above definitions to 



A(q) = 7 + A 



A x(q) 

1 - Aox(q) 



(24a) 



We see that the high-dimensional irreducible vertex is 
completely determined from a single local (mean-field) 
parameter Ao and the two-particle bubble %(q). Sum- 
ming both sides of Eq. I|24a|l over momenta we obtain an 
equation for the local two-particle irreducible vertex 



An 



A 



x(q) 



-y- 

N^l- Aox(q) 



(24b) 



Knowing the local part of the two-particle irreducible ver- 
tex Ao we can reconstruct the full two-particle vertex in 
high spatial dimensions. We have 



kk 



(q) = 7 



Ar 



Aox(q) 



A x(k + k' + q) 



_ 1 - Aox(q) 1 - A x(k + k' + q) 



(25) 
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where, in analogy to the non-self-consistent high dimen- 
sional vertex, Eq. JTHJ, A = A /(l + A G + G_) and 
X(q)=x(q)-G + G_. 

Equations (|23[1 - l|25[l form an approximation for two- 
particle functions with a single, local parameter deter- 
mined self-consistently. Such an approximation is a self- 
consistent extension of the high-dimensional limit of the 
two-particle vertex, Eq. (|18|) . We hence can call it a 
mean-field approximation for two-particle functions of 
noninteracting disordered electron systems. The self- 
consistently determined mean-field parameter Ao exactly 
reproduces the leading 1/d correction to the CPA irre- 
ducible vertex A. The fundamental self-consistent equa- 
tion of this approximation, Eq. (|24b() . has a typical mean- 
field character. That is, it can be used in any dimension. 
The lattice structure enters the equation only through 
the momentum summation running over the first Bril- 
louin zone. Notice, however, that we cannot reduce 
the momentum summation in two-particle functions to 
an integral over the density of states. In high but fi- 
nite spatial dimensions we have to sum over momenta in 
Eq. (|24b(l by using the Gaussian rules only asymptoti- 
cally for d — > oo 

D. Self-energy and diffusion pole in high 
dimensions 

The mean-field theory constructed in the preceding 
subsection is a self-consistent approximation for the two- 
particle vertex where one-particle propagators are as- 
sumed to be external functions. Such a situation cannot 
be the final stage of the theory, since due to conservation 
laws and thermodynamic consistence the one- and two- 
particle functions are correlated. Actually, the electron- 
hole irreducible vertex is connected with the self-energy 
for noninteracting disordered electrons via the Vollhardt- 
Wolfle Ward identityii 

E*(k, E + u) - £*(k, E) = I A ^( E + ^ E ) 

k' 

x [G R {\s!,E + w) -G A {k',E)] . (26) 

Here we denoted Y, R , E A the retarded and advanced self- 
energy and Af$(E+w,E) = A^,(E+uj+iO + ,E-iO + ;0). 
Note that the irreducible vertex A in Eq. (|26|) does not 
have bar, i. e., it is defined via the Bethe-Salpeter equa- 
tion with the full one-electron propagators G R , G A . The 
Ward identity says that if we modify the electron-hole 
irreducible vertex we have to change adequately the one- 
electron self-energy and vice versa. We hence cannot 
approximate independently the irreducible vertex with- 
out changing appropriately the self-energy. If we have an 
analytic expression for the self-energy as a functional of 
the one-electron resolvent, we can use a differential form 
of the Ward identity to determine the irreducible 
vertex. Then the vertex function is determined directly 
from the one-electron propagator. Such a construction of 



the irreducible vertex does not lead to bifurcation points 
and multiple solutions, i. e., to a phase transition, un- 
less we find them in the self-energy itself. It is normally 
very difficult to determine bifurcation points in the self- 
energy that is a bounded function. It is more convenient 
to search for possible bifurcation points in two-particle 
functions that may display divergences. 

It is clear that we need self-consistent equations for 
functions that could have multiple solutions. Such a self- 
consistent approximation for the irreducible vertex was 
achieved in the preceding subsection. A self-consistent 
equation for the irreducible vertex cannot be derived 
from a self-energy directly. To achieve consistence be- 
tween the one- and two-particle functions in this case we 
have to determine the self-energy as a functional of the 
irreducible vertex. This must be done in concord with 
the Ward identity 

The self-energy as a functional of the irreducible ver- 
tex is overdetermined from identity l|26|) . We can nev- 
ertheless use the Vollhardt-Wolfle identity to determine 
the self-energy from the vertex function as suggested in 
Ref. |52 and used in the parquet approach from Ref. 
We use only a specific element of the Vollhardt-Wolfle 
identity to determine the imaginary part of the self- 
energy. In the mean-field case, where the irreducible ver- 
tex is local, we can write a generalized CPA equation 

%Z R (E) = A RA (E,E)%G R (E) . (27a) 

Consistence, or negative definiteness of the self-energy, 
demands that the local element of the irreducible vertex 
A RA (E, E) = N- 2 £ kk , A kk ,(E + iO+ E - i0+; 0), de- 
termined from the mean- field equation Ij24bf) . is positive. 

We cannot find the real part of the self-energy directly 
from the irreducible vertex. Instead, we use causality 
of the self-energy and the Kramers-Kronig relation ex- 
pressing the real part of an analytic function as a Hilbert 
transform of its imaginary part. We have 

mE)=Sx+P r^^m. (27b) 

J-oo I" E'-E 

Equations l)27|) now determine the self-energy from the 
mean- field (local) irreducible vertex Ao. Notice, how- 
ever, that to determine the self-energy at one frequency 
we have to know the irreducible vertex in the whole fre- 
quency range. Equations i|27|) complete the mean-field 
theory for vertex functions, Eqs. 1)23(1 and Eqs. Q24J1 . and 
make it a consistent approximation with proper analytic 
properties of one- and two-particle Green functions. 

The Vollhardt-Wolfle identity, its specific form from 
Eq. H27a|) . not only serves as a means for a consistent 
determination of a causal self-energy from a given irre- 
ducible vertex, but it is also indispensable for the exis- 
tence of the diffusion pole in the electron-hole correla- 
tion function. We now show in what form the diffusion 
pole survives in the mean-field theory for vertex functions 
with the self-energy determined by Eqs. I|27l) . 
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The electron-hole correlation function is defined from 
the averaged two-particle Green function as 



1 



*r(q^) = ^^4k'(£+^+*0 + ,£-*0+;q) . (28) 

kk' 

The two-particle Green function is evaluated with the full 
two-particle vertex via Eq. Q . 

From Eq. l(2Ta|) and from x flj4 (0) = QfG^/SSE* we ob- 
tain that both the denominators in the high-dimensional 
limit of the two-particle vertex T RA vanish at zero trans- 
fer momenta. Hence, the mean-field approximation 
for vertex functions contains the diffusion pole in the 
electron-hole channel (first term in the brackets on the 
r.h.s. of Eq. an d the Cooper pole in the electron- 

electron channel (second term). Only the diffusion pole 
survives as a singularity in the electron-hole correlation 
function, Eq. (J28J. 

To derive the low-energy behavior of the electron-hole 
correlation function we denote 



A, 



w=0 



and 



D%{uj) = 25iT± u (E)A$ A (E+u>,E) 



d X HA (<l) 



d(q 2 



(29a) 



(29b) 



g=0 



With these two definitions we find the high-dimensional 
asymptotics of the low-energy limit of the electron-hole 
correlation function at zero temperature to be 



1>^(q,u,) 



2nni 



-iA E u + D E {u)c? 



(30) 



where tie is the density of states at the Fermi energy E. 

The low-energy asymptotics of the electron-hole cor- 
relation function serves as an important tool for testing 
consistence of approximations. We find from gauge in- 
variance and the (unrestricted) Ward identity (|26|l that 
the electron-hole correlation function should display the 
diffusion pole in form of Eq. I|3UI) with A% = 1 for ar- 
bitrary disorder strength. It then means that the low- 
energy behavior of the electron-hole correlation func- 
tion is controlled by a single parameter, the bare dy- 
namical diffusion constant D e (uj)H^ However, we al- 
ready found in Ref. ^] that the constant A e increases 
with the disorder strength and becomes infinite at the 
Anderson localization transition. The disorder depen- 
dent weight of the diffusion pole ue/Ae is the cen- 
tral unexpected feature of the mean-field theory for ver- 
tex functions. This mean- field approximation obeys the 
Vollhardt-Wolfle identity only in the limit to zero fre- 
quency, Eq. (|27a|l . and not for finite energy differences. 
We found a consistent explanation for the decrease of 
the weight of the diffusion pole with increasing disorder 
strength^ The weight of the diffusion pole, tie/Ae, ex- 
presses a portion of extended (diffusive) states from all 
available states at the Fermi energy determined by the 



density of states calculated from the one-electron Green 
function, n E = -%G r (E)/t: . 

The dependence of the weight of the diffusion pole 
on the disorder strength found in the mean-field theory 
for averaged two-particle functions could be an artifact 
of this specific approximation. We could still hope that 
the full exact solution recovers the invariant weight of 
the diffusion pole expected from the unrestricted con- 
servation laws for averaged Green functions. Based 
only on approximation-free arguments we found that the 
Vollhardt-Wolfle identity (|26() for finite frequencies can- 
not be fulfilled in any finite dimension if the electron-hole 
irreducible vertex A eh contains the Cooper pole. Enforc- 
ing the full form of the Vollhardt-Wolfle identity with 
the Cooper pole in the electron-hole irreducible vertex 
inevitably leads to a self-energy being a nonanalytic func- 
tion of its energy argument for almost all Fermi energies 
within the energy bands^iSi It hence seems that the 
high-dimensional behavior of the two-particle vertex and 
the disorder-dependent weight of diffusion pole, disclosed 
by the asymptotic mean-field solution, are generic fea- 
tures of the Anderson model of disordered electrons. At 
least for theories that can be analytically continued from 
the limit to high spatial dimensions. 



IV. DISCUSSION AND CONCLUSIONS 

We presented in this paper two ways how to reach 
mean-field-like approximations for noninteracting disor- 
dered electron systems. The first one, being the standard 
thermodynamic mean-field theory known from many- 
body systems, uses the limit to infinite spatial dimen- 
sions applied to the generating, configurationally aver- 
aged thermodynamic potential. The limit to infinite spa- 
tial dimensions enables one to separate the diagonal and 
off-diagonal elements of the one-particle propagator and 
its self-energy and to find an explicit representation for 
the generating functional. The local one-particle func- 
tions from the generating functional serve as generalized 
variational parameters, the physical values of which are 
determined from stationarity equations for the generat- 
ing functional. The higher-order Green and vertex func- 
tions are determined from responses of the system to local 
external perturbations. In this way the construction of 
a mean-field approximation is consistent and unambigu- 
ous in the determination of one-particle functions as they 
are the only ones entering the generating functional. The 
higher-order Green functions are defined uniquely only in 
their local parts. 

The nonlocal parts of two-particle Green functions are 
no longer determined from the local thermodynamic the- 
ory uniquely. We can either use the standard construc- 
tion of Baym and Kadanoff or we can directly apply the 
asymptotic limit to infinite spatial dimensions to two- 
particle functions. The outcome of these two construc- 
tions is not identical. In the former way we miss some 
of the leading-order (maximally crossed) diagrams and 
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lose the electron-hole symmetry at the two-particle level. 
These deficiencies severely discredit reliability of the ther- 
modynamic mean-field theory in the calculation of spa- 
tial coherence and transport properties of disordered sys- 
tems. To overcome these drawbacks we proposed another 
route toward a mean-ficld-like approximation for vertex 
functions based on a direct analysis of two-particle func- 
tions in high spatial dimensions. 

The incapability of the thermodynamic mean-field 
theory to correctly describe nonlocal correlations in two- 
particle functions is caused by the degeneracy of local 
theories with elastic scatterings. Static local approxi- 
mations are unable to distinguish between electrons and 
holes. Only quantum dynamics or multiple scatterings 
on spatially distinct impurities can discern the motion of 
an electron from the motion of a hole. The distinguisha- 
bility of electrons and holes is of principal importance for 
encountering backscattering effects and for a two-particle 
self-consistence used in the construction of a mean-field 
theory for vertex functions. 

A mean-field approximation for vertex functions in 
disordered electron systems was constructed from the 
asymptotic limit to high spatial dimensions, where, unlike 
the thermodynamic mean field, the lattice dimension is 
high but finite. Alike the strict limit d = oo, the asymp- 
totic behavior in high dimensions leads to significant sim- 
plifications in momentum convolutions that enable us to 
reduce the approximation to a mean-ficld-like one with 
a local generator determined from a self-consistent equa- 
tion. We applied the asymptotic limit directly within a 
diagrammatic expansion around the CPA for two-particle 
functions. The basic ingredients for the construction 
were parquet equations for the two-particle irreducible 
vertices from the electron-hole and the electron-electron 
scattering channels. Using the electron- hole symmetry 
at the two-particle level and the asymptotic limit to high 
dimensions we succeeded in producing a self-consistent 
approximation for the local part of the electron-hole ir- 
reducible vertex. It is a self-consistent 1/d extension of 
the CPA irreducible vertex. The mean-field theory for 
vertex functions then determines in a unique way the full 
two-particle vertex that correctly reproduces the limit to 
infinite dimensions with the electron-hole symmetry at 
both one- and two-particle levels. 

The mean-field theory for vertex functions is an ap- 
proximation for two-particle functions. The one-electron 
functions, used as an input for the two-particle equations, 
are then calculated from the vertex function via a spe- 
cific form of the Vollhardt-Wolfle Ward identity and the 
Kramers-Kronig relation. With this extension of the the- 
ory to one-particle functions we accomplished an approx- 
imate description of disordered systems fulfilling all con- 
sistency conditions on one- and two-particle functions. 

The most important achievement of this mean-field 
theory is its ability to describe the disorder-induced van- 
ishing of diffusion, that is, the Anderson localization 
transition. This theory succeeded for the first time to 
bridge qualitatively correctly the weak and the strong 



disorder limits and to cover the split-band limit as well 
as vanishing of diffusion. The other existing approaches 
have concentrated on only one of the two phenomena. 
They either miss the two-particle self-consistence or do 
not consistently match the one- and two-particle func- 
tions. 

The consistence between the one-electron self-energy 
and the electron-hole irreducible vertex is essential for 
credibility of approximate treatments of the Anderson 
metal-insulator transition. Only with this relation cor- 
rectly taken into account we obtain the proper form of 
the diffusion pole and electron diffusion on long distances. 
In this respect the mean-field theory for vertex functions 
leads, surprisingly against the common expectations, to 
a nonconserving weight of the diffusion pole and its de- 
pendence on the disorder strength. The thermodynamic 
mean-field theory and all other approaches to the Ander- 
son localization transition assume or use the unrestricted 
form of the Ward identity (|26|l being a consequence of 
conservation of the norm of the wave function for all 
configurations of the random potential. We argued al- 
ready earlier that the Hilbert space of Bloch waves is 
incomplete in the sense that it does not encompass the 
eigenstates of all configurations of the random potential. 
At a given energy we observe even in the metallic regime 
macroscopically relevant numbers of configurations with 
localized states^ 

Vanishing of the diffusion pole at the Anderson metal- 
insulator transition and in the localized phase modifies 
the existing picture of the critical behavior for vanishing 
of diffusion. In the standard approaches, such as the non- 
linear sigma model or the Vollhardt-Wolfle self-consistent 
approximation, the dynamical diffusion constant is the 
only parameter controlling the low-energy behavior of 
the electron-hole correlation function. In our mean-field 
theory we have apart from the diffusion constant also 
the weight of the diffusion pole that significantly influ- 
ences the description of the long-range correlations and 
diffusion. In the critical region, however, only the con- 
stant Ae — > oo from Eq. (|29a|l is relevant and all critical 
scales can be derived from it. For instance the renor- 
malized diffusion constant De = D e /Ae vanishes at the 
localization transition with the diverging parameter Ae, 
etc. It means that the Anderson localization transition 
is compatible with a one-parametric scaling theory. 

Although the one-parameter scaling holds for the An- 
derson localization transition, two relevant parameters 
in the critical region, ue/Ae and D e /Ae, nevertheless 
lead to a modification of the critical behavior deduced 
from the field-theoretic approaches. The two parame- 
ters stand for two quantities influencing the electrical 
conductivity. The former expresses the number of ex- 
tended states at the Fermi energy E and the latter the 
averaged velocity of the diffusive particles. Both quan- 
tities go simultaneously to zero at the Anderson metal- 
insulator transition. It is straightforward to verify that 
the mean-field theory for vertex functions predicts that 
Ae ~ |A C — A| -1 / 2 , where A is the bare disorder strength 
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and A c its critical value. Having two vanishing parame- 
ters we have to distinguish two types of the critical be- 
havior. First, we have properties of individual electrons. 
One of them is diffusion as seen from the semiclassical dif- 
fusion equation. This is quantitatively described by the 
renormalized diffusion constant De ~ |A C — Al 1 / 2 . Sec- 
ond, we have statistical values describing the disordered 
sample as a whole. Among them the averaged conduc- 
tivity is the most interesting one. It is proportional to 
the product of the number of available diffusive states 
and the renormalized diffusion constant, so that we have 
a ~ ueD e /A 2 ~ | A c — Al 1 . Notice, however, that there 
is not a direct relation between the diffusion constant 
and the conductivity calculated from the Kubo formula, 
since due to deviations from the Ward identity, the Ein- 
stein relation does not hold. Finally we can also deduce 
the critical exponent for the localization length in the lo- 
calized phase. Its square is inversely proportional to the 
order parameter, being the imaginary part of the local 
irreducible vertex QA RA (E + 0+, E - 0+) - |A - \ c \ 1/2 . 
We hence obtain £ ~ |A — A c | -1 ^ 4 . The critical expo- 
nent for the conductivity equals the mean-field exponent 
from the Vollhardt- Wolfle theory but the critical expo- 
nent for the localization length is one half of their value. 
We remind that the critical behavior obtained from the 
mean-field theory for vertex functions holds only for di- 
mensions d > d u = 4. We expect corrections to the 
mean-field critical behavior in lower dimensions. 

To conclude, we demonstrated that to reach a reli- 
able quantitative description of the Anderson localiza- 
tion transition one has to sum up self-consistently non- 



local contributions to the electron-hole and the electron- 
electron irreducible vertex functions. It can be achieved 
in a mean-field manner via the asymptotic limit to high 
spatial dimensions leading to a self-consistent extension 
of the CPA local vertex. The most important conclusion 
of this mean-field theory of Anderson localization is that 
the weight of the diffusion pole is not conserved and that 
the diffusion pole is absent in the localized phase. The 
decreasing weight of the diffusion pole with the increas- 
ing disorder strength is a consequence of incompleteness 
of the Hilbert space of Bloch waves. At any Fermi energy 
there are macroscopically relevant numbers of configura- 
tions with localized as well as with delocalized states. 
The number of configurations with extended states de- 
creases with increasing the disorder strength and vanishes 
at the localization transition. This feature can essentially 
be checked by other, e. g. numerical, means. Due to the 
disorder-dependent weight of the diffusion pole the Ein- 
stein relation does not hold and we have to distinguish 
individual and statistical transport properties of the dis- 
ordered systems. One should have this in mind when 
calculating the critical behavior of the Anderson local- 
ization transition. 
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